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S" 

A trial wave function for two-dimensional quantum dot helium in an arbitrary perpendicular 
magnetic field (a system of two interacting electrons in a two-dimensional parabolic confinement 
£S) . potential) is introduced. A key ingredient of this trial wave function is a Jastrow pair correlation 

factor that has a displaced Gaussian form. The above choice of the pair correlation factor is instru- 
mental on assuring the overall quality of the wave function at all values of the magnetic field. Exact 
numerical diagonalization results are used to gauge the quality of the proposed trial wave function. 
We find out that this trial wave function is an excellent representation of the true ground state at 
(■ — ' all values of the magnetic field including weak (or zero) and strong magnetic fields. 

PACS numbers: 73.21. La, 31.15.Pf 
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INTRODUCTION 



Among the variety of two-dimensional (2D) few electron quantum dots [H, 0, S, S, i, S 0, B S E3, EH, E3, E, Q , the 
two-electron (N — 2) quantum dot stands out as a truly remarkable system that, despite the simplicity, shows very 
rich phenomena and posseses characteristic features which persist to larger systems. Such a system oftenly refered 
to as 2D quantum dot helium exhibits a highly complex behavior in presence of a perpendicular magnetic field. Its 
ground state energetics is unusually complicated and intricate singlet to triplet spin state transitions occur as the 
I , magnetic field is varied. [l5|, 

The Hamiltonian of 2D quantum dot helium in a perpendicular magnetic field can be written as: 
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■ where p i = (fiixTPiy) and Pi — (pix,Piy) are, respectively, the 2D momentum and position of the i-th electron, m is 
electron's mass, — e(e > 0) is electron's charge, g e is electron's g-factor, p,B is Bohr's magneton, e r is the dielectric 
00 ' constant, Huq is the parabolic confinement energy, S z is the z-component of the total spin, B z is the perpendicular 
magnetic field, Li Z is the z-component angular momentum operator for the z-th electron and uj c = eB z /m > is the 
cyclotron frequency. 

There have been several studies of 2D quantum dot helium with or without magnetic field employing a variety 



of techniques. For example, Merkt et al. [17[ studied the energy spectra of two interacting electrons in a parabolic 



potential in absence and in presence of a perpendicular magnetic field employing the exact numerical diagonalization 



technique. Wagner et al. [15| used the exact numerical diagonalization technique to study 2D quantum dot helium 
in a perpendicular magnetic field and to predict oscillations between spin-singlet and spin-triplet states as a function 
of the magnetic field strength. Pfannkuche et al. compared energies, pair correlation functions, and particle 
densities obtained from the Hartree, Hartree-Fock (HF) , and exact numerical diagonalization method and pointed out 
the unsuitability of the HF approach at weak magnetic fields. Pfannkuche et al. used the diagonalization data 
to formulate a theory of 2D quantum dot helium in a perpendicular magnetic field, explaining how the competition 
between the Coulomb interaction and the binding forces due to confinement and the magnetic field induces ground 
state transitions. Harju et al. [2(| introduced a recipe of how to build a trial wave function for 2D quantum dot 
helium in a perpendicular magnetic field and applied the variational Monte Carlo (VMC) technique. His Jastrow pair 
correlation factor has two variational parameters and takes into account the mixing of different Landau levels for the 
relative motion, though the ansatz does not have the right asymptotic behavior at large interparticle distances. 
A wealth of information in 2D quantum dots, in general, and 2D quantum dot helium in particular, is provided 



by the rotating Wigner (or electron) molecule (RWM or REM) theory. [2l|, |22| Specific REM trial wave functions 
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have been recently derived [23[ for 2D quantum dot systems in high magnetic fields. The REM wave functions are 
constructed by first breaking the rotational symmetry at the unrestricted Hartree-Fock (UHF) level and, secondly, 
restoring the circular symmetry via post HF methods and projection techniques. The UHF-level REM wave 
function describes a Wigner molecule that is considered as a rigid rotor, while the second step of restoring the circular 
symmetry implies rotations of such molecules. At zero and weak magnetic fields the broken-symmetry UHF orbitals 
need be determined numerically while at high magnetic fields they can be well approximated by one-particle lowest 
Landau level (LLL) Gaussian functions that are centered at positions that correspond to the classical equilibrium 



configuration of N point charges in a harmonic trap. Yannouleas and Landman [22j used such RWM trial wave 
function to study 2D quantum dot helium in a magnetic field. They calculated the ground state energies as a function 
of the magnetic field for the spin triplet states of 2D quantum dot helium and noted that for high magnetic fields ( 
larger that 7 T, see Fig. 1 in Ref. [22j) the RWM and exact diagonalization results are practically the same. 

Some of the main theoretical methods used to study quantum dots are: analytical calculations [24l . I251 1261 . 127 



exact numerical diagonalizations 



28, ]2J,j30|, l3l(, quantum Monte Carlo (QMC) methods 



32, |33|, 13J, 123], density 
44] These methods have 



functional theory methods 2fl 22, 23, El, So] and Hartree-Fock mean- field theory. (U 42, 43 
advantages and disadvantages, most notably some methods that deliver good results for the ground state properties 
may be inadequate in describing excited state properties and vice versa. Here we are limiting our discussion to the 
calculation of ground state properties only. Considering the limitations of several methods on obtaining ground state 
properties (for example: exact numerical diagonalization methods are applicable only for few electrons, Hartree-Fock 
and perturbation theory lack the desired accuracy, etc), the use of QMC methods, such as the VMC technique seems 
to be the best strategy in the long run. Therefore the quest for better, yet simple, trial ground state wave functions is 
always timely. A high quality trial ground state wave function is essential not only for the VMC method, but also for 
the more sophisticated diffusion Monte Carlo (DMC) method that relies on a guiding trial wave function. Compared 
to other methods, QMC methods have the greatest advantage of all, since they can be extended to larger number of 
electrons in a straightforward manner and are very accurate. 

In this work we introduce a novel trial wave function to describe the ground state of 2D quantum dot helium in a 
perpendicular magnetic field. Such a wave function is written as a product of a Laughlin-type wave function 4^, 46|, 47] 
with a Jastrow pair correlation factor and has the form: 



*G5i,/&) = J(p ia ) x ( Zl - z 2 )l^l exp (-^2^) . (2) 
where the Jastrow factor, J(pi2) is: 

J(pi 2 ) = exp f-yPi2 +c 6/3i2 V (3) 

The 2D position coordinate, Zj — Xj —iyj of the Laughlin component of the wave function is given in complex notation, 
P12 = \Pi — /0*2 1 is the inter-electron distance, m z = \m z \ = 0, 1, . . . is the angular momentum number, b and c are 
non-negative variational parameters to be optimized. We have In = y f fi/(2mQ) and SI 2 = luq + {lo c /2) 2 . The effective 



magnetic length, Iq reduces to the electronic magnetic length, la = y ^g- when there is no confinement (ojq = 0) or 

when the magnetic field is very large ( ui c /uio — > oc ). The effective magnetic length, Iq can be written in terms of 
the inverse oscillator length, a — y^muo/h as: 




(4) 

In this way it is easy to recover the 2D harmonic oscillator states in the limit of zero magnetic field ( lo c /ljq — > 0). 
The parity of the ground state space wave function depends on the value of angular momentum \m z \. For even 
values, m z \ — 0, 2, 4, . . . the space wave function is symmetric and the spin function corresponds to a spin-singlet 
state (S = 0), while for odd values, \m z \ = 1, 3, 5, . . . the space wave function is antisymmetric and the spin function 
becomes a spin-triplet (S — 1). 

What makes this trial wave function rather unique is its dis plac ed Gaussian Jastrow pair correlation factor, J{p\2) 1 
which is different from earlier choices in the literature. 2(J 33, 48| The rationale behind the choice displaced Gaussian 



choice of correlation factor can be better understood if one considers a pair of electrons in zero magnetic field. In 
absence of an electronic repulsion between the electrons the relative coordinate ground state wave function with be 
a Gaussian centered at coordinate, p\2 = 0, will have zero angular momentum and will correspond to a spin singlet 
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state. With the Coulomb repulsion the ground state will still have zero angular momentum [33J therefore it is plausible 
to expect that the main effect of the Coulomb correlation is simply to further separate the electrons resulting in a 
new relative coordinate ground state wave function which will resemble a Gaussian centered at P12 7^ values. The 
choice in Eq.Q mimicks this physical effect. 

Trial wave functions, such as the REM wave functions (for N = 2 and N > 2 electrons) also take into consideration 
the relative separation of electrons and have led to a rich physics regarding Wigner crystallization and the rotation 
of the electron molecules formed in high magnetic fields. [22|, |49( Within the framework of the REM wave functions, 
which by construction are crystalline in character, the separation between electrons is achieved at the one-particle 
level right at the start. For instance, in the high magnetic field regime, the REM wave function's Slater determinants 
contains LLL one-particle Gaussian orbitals which are centered at different positions, Zj (in complex notation) and 
have the form: 



u(z,Zj) = exp 



4/ 2 



(5) 



where Zj coincide with the classical equilibrium positions of classical point charges and Iq is the electronic magnetic 
length. In the case of our variational wave function, the separation between electrons is achieved at the two-particle 
(pair) level through the displaced Gaussian pair correlation factor. Despite the same idea of optimizing the separation 
between electrons, the displaced Gaussian pair correlation factors considered in this work have a two-body structure 
which makes it different from the one-particle (displaced) Gaussian functions that approximate the one-particle UHF 
orbitals of the REM wave functions at high magnetic fields. [23| 

Despite the usual controversies involved on the selection of a Jastrow correlation factor it appears clearly that the 
displaced Gaussian pair correlation factor has all the attributes to capture effectively most of the electronic correlations 
present on the ground state of such system. Some indication of the eventual high quality of the trial wave function 
also comes from a previous study of 2D quantum dot helium at zero magnetic field [501 ] where the same Jastrow pair 
correlation factor was used. 

Motivated by these arguments, it is the objective of this work to perform a complete study of 2D quantum dot 
helium in an arbitrary perpendicular magnetic field that ranges from very weak (or zero) to very strong using the trial 
wave function introduced in Eq.Q and Eq.([2]). In the process we will test how effective the displaced Gaussian pair 
correlation factor is on representing the effect of electronic correlations at arbitrary magnetic fields. Exact numerical 
diagonalization calculations will be used as a gauge of accuracy of the variational results. 

The paper is organized as follows: In Section we introduce the variational method, the trial wave function and 
calculate various quantities corresponding to such a wave function. In Section we show the results obtained from the 
variational wave function and compare them to the exact diagonalization method. A discussion of the results is given 
in Section and concluding remarks can be found in Section . 



VARIATIONAL THEORY 



In this section we apply the variational method to study 2D quantum dot helium in a perpendicular magnetic field 
using the trial wave function of Eq.([2]) with the Jastrow pair correlation factor having the displaced Gaussian form 
given in Eq.Q. We will show that, after optimization, the proposed trial wave function is an excellent representation 
of the true ground state at any value of the magnetic field and compares very favorably to the exact numerical 
diagonalization results. 

There are two dimensionless parameters that determine the behavior and properties of the system under consider- 
ation: 



47reoe r ' loq ' 

where A gauges the strength of the Coulomb correlation relative to the confinement energy and 7 measures the strength 
of the magnetic field relative to confinement. One cam immediately see that A = I /as, where I = 1/a is the harmonic 
oscillator length and as = 47reoe r 7i 2 / (me 2 ) is the effective Bohr radius. 

Since the parity of space wave function is determined by the value of the angular momentum, \m z \, the ground 
state angular momentum value determines whether the ground state corresponds to a spin-singlet or spin-triplet state. 
Therefore, a stringent test of quality for this trial wave function is to check whether the lowest variational energy state 
has always the correct angular momentum number as calculated from the numerical diagonalizations under different 
combinations of Coulomb correlation, confinement and magnetic field. 
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In a general situation where both Coulomb correlation and magnetic field are present, there is no way to anticipate 
the correct value of the ground state angular momentum number. An exception are the simple cases of: (i) absence 
of Coulomb correlations or (ii) absence of a perpendicular magnetic field, where it is straightforward to prove that 
the ground state is expected to have zero angular momentum. 

Obviously, such behavior should be reflected by the trial wave function under investigation. With no Coulomb 
correlation (A = 0), the ground state has zero angular momemtum (\m z \ — 0) both in presence or absence of the 
perpendicular magnetic field. Under these conditions, the corresponding trial wave function becomes that of Eq.([2]) 
with Jastrow correlation J(p±2) = 1 and angular momentum, \m z \ — as expected. In absence of a perpendicular 
magnetic field, the groundstatc still has zero angular momentum, (\m z \ = 0) with or without the Coulomb correlation, 
therefore a trial wave function with J{pri) ^ 1 and angular momentum, \m z \ — is again consistent with the expected 
scenario. 

However, when both Coulomb correlation and magnetic field are present, the situation changes drastically. As the 
magnetic field increases, a groundstate with nonzero angular momentums (\m z \ ^ 0) may arise, therefore in addition 
to the Jastrow pair correlation factor also the Laughlin factor starts contributing on keeping the electrons apart in a 
more effective way. 

A calculation of the expectation value of the Hamiltonian with respect to the trial wave function, E — 
gives: 



E 

e(B,c,i,m z ,X) = = — 



\m z \ 



■7 + 



-B 2 h(B, c, 7 , m z ) + (1 + £ ) MI l^r ) + XB fa(B, c, 7, m t ) 



f(B,c,~/,m z ) 



7 2 

1 + T 



(7) 



where B = b/a and c are variational parameters, 7 = lo c /ujq is linearly proportional to the magnetic field (oc B z ) 
and A = e 2 a / 1 (Aireae r fujjQ) is the Coulomb correlation parameter. All the above quantities are dimensionless and the 
ground state energy is given in units of Tiujq. Note that the Zeeman energy term is not specifically included in the 
expression for the variational energy. 

The functions: /i,2,3 and /, as well as the function g, depend on the variables specified on their arguments and in 
integral form are given by: 



/i(B,c,7,m z ) = r o °° dttg(t,B,c^,m z )y. 



til 



1 

2W 1 



1 + 



Kl 1 



f 2 (B,c,j,m z ) = r o °° dtt 3 g(t,B,c,j,m z ) , 
f 3 (B,c,j,m z ) = J °° dtg(t,B,c,j,m z ) , 
f(B,c,j,m z ) = r o °° dttg(t,B,c,j,m z ) , 



3 (i,B,c,7,m z ) = t 2 l^l 



e 2B- 



-r+2ct 



1+^-2 



(8) 



where t = b P12 is an auxiliary variable introduced to simplify the calculation of integrals. The optimization procedure 
is straightforward: given the values of Coulomb and magnetic field parameters, A and 7 we calculate the lowest 
energies for a set of integer values of m z by optimizing the variational parameters, B, and c through standard 
numerical procedures. 



RESULTS 



The best way to gauge the accuracy of the trial wave function is to directly compare the variational results to 
exact numerical diagonalization values. Table U displays the diagonalization ground state energies, e = E/(JujJq) for 
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TABLE I: The exact numerical diagonalization ground state energies, e = E/(hu)o) for 2D quantum dot helium subject to 
a perpendicular magnetic field as a function of dimensionless Coulomb coupling parameter, A — e 2 a/(Aneo^rfi^ l o) — 0, 1, ... 6 
and values of magnetic field, j = uj c /uo = 0, 1, . . . , 5. The angular momentum, m z of the ground state is also specified. The 
parameter a = ^Jmuo/h has the dimensionality of an inverse length. 





7 = 


7 = 1 


7 = 2 


7 = 3 


7 = 4 


7 = 5 


A = 


2.00000 

m z = 


2.23607 

m z = 


2.82843 
m z = 


3.60555 

m z = 


4.47214 

m z = 


5.38516 

m z = 


A = 1 


3.00097 
m z — 


3.30508 

m z = 


3.95732 
m z = 1 


4.71894 

m z = 1 


5.61430 
m z = 1 


6.53067 

m z = 2 


A = 2 


3.72143 

m z = 


4.06684 
m z = l 


4.61879 

m z = 1 


5.43123 

m z = 2 


6.30766 

m z = 2 


7.22681 

m z = 3 


A = 3 


4.31872 

m z = 


4.60594 
m z = 1 


5.23689 
m z = 1 


6.01256 

m z = 2 


6.89002 

m z = 3 


7.81384 

m z = 4 


A = 4 


4.84780 
m z = 


5.11165 

m z = 1 


5.73642 

m z = 2 


6.53522 

m z = 3 


7.41600 

m z = 4 


8.33874 

m z = 5 


A = 5 


5.33224 
m z — 


5.58995 
m z = 1 


6.21499 

m z = 2 


7.01716 

m z = 3 


7.90109 

m z = 4 


8.82281 

m z = 6 


A = 6 


5.78429 
m z = 


6.04534 

m z = 1 


6.67999 

m z = 2 


7.46782 

m z = 4 


8.34530 

m z = 5 


9.27057 

m z = 6 



2D quantum dot helium in a perpendicular magnetic field for values of Coulomb correlation, A = 1, ... 6, and values 
of magnetic field, 7 = 0, 1, . . . , 5. The ground state angular momentum, m z is also specified. 

In Table HT1 we show the variational ground state energies and optimal values of parameters, B and c for 2D quantum 
dot helium in a perpendicular magnetic field at different A-s and 7-s. The results are rounded in the last digit. The 
variational energies shown in Table [TT] are in excellent agreement with numerical diagonalization results reported 
in Table |U This agreement holds for the whole range of Coulomb correlations and perpendicular magnetic fields 
considered in this work. 

In the strong magnetic field limit, the variational energies are practically identical (within the range of very small 
statistical errors) to the exact numerical diagonalization values indicating the overall excellent quality of the trial wave 
function. Even more remarkable is the fact that for any combination of A-s and 7-s the angular momentum of the 
lowest variational energy state is always reached at the exact value obtained from the exact numerical diagonalizations. 
We note there are cases, such as A = 5 and 7 = 4 where the energy difference between ground state and higher states 
with different angular momentum is extremely small. ( For A = 5 and 7 = 4 the diagonalization method gives a 
ground state with energy: e = 7.90109 and angular momentum \m z \ = 4, while the first excited state has an energy 
slightly higher: e = 7.90640 and angular momentum \m z \ = 5 ). Nevertheless, in all occasions the trial wave function 
with lowest energy has an angular momentum that corresponds to the exact diagonalization value. 

Both diagonalization and variational results are in full agreement and confirm the expected outcome that: (i) in 
absence of Coulomb correlations (A = 0) or (ii) in absence of a perpendicular magnetic field (7 = 0), the ground 
state has zero angular momentum. However, when both A / and 7 7^ 0, ground states with nonzero angular 
momentum arise. For very large values of the perpendicular magnetic field the ground state has increasingly large 
angular momentum values where each change of |m z | indicates a singlet to triplet spin state transition, a phenomena 
that has been observed in recent experiments [16j with ultrasmall quantum dots. Since the ground state spin of 2D 
quantum dot helium can be either singlet (S — 0) or triplet (S — 1) it is plausible to expect that this quantum dot 
has the potential to serve as qubit of a quantum computer, with the magnetic field tuning the transition between the 
two spin states, an idea suggested by Burkard et al. [51| 

Because the change of angular momentum indicates a spin state transition, it is the mean square distance between 
the two electrons (which is directly related to the angular momentum number) that should indicate jumps or non- 
monotonic behavior contrary to the dependence of energy on A that is monotonic. Therefore, in addition to ground 
state energies, we also calculated the mean square distance between two electrons, (\pi — P2I 2 ) for a wide range of 
A-s and 7-s. The results are displayed in Table. IIIII where we show a 2 (\pi — P2I 2 ) for values of A = 0, 1, . . . , 10 and 
7 = 0, 1, . . .5. In absence of Coulomb correlations, the increase of the magnetic field brings electrons closer to each 
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TABLE II: The variational ground state energies, e = E/(huo), angular momentum values, m z , as well as optimal parameter 
values, B = b/a and c for 2D quantum dot helium subject to a perpendicular magnetic field as a function of dimensionless 
Coulomb coupling parameter, A = e 2 a/(47reo£r?ii^o) = 0, 1, ... 6 and values of magnetic field, 7 = lo c /^o = 0, 1, . . . , 5. The 
parameter a = t/muj/S has the dimensionality of an inverse length. 





7 = 


7 = 1 


7 = 2 


7 = 3 


7 = 4 


7 = 5 


A = 





2.00000 


2.23607 


2.82843 


3.60555 


4.47214 


5.38516 


m z 






















B 






















c 






















A = 


1 


3.00174 


3.30578 


3.95737 


4.71899 


5.61435 


6.53068 


m z 










1 


1 


1 


2 


B 




0.40185 


0.41627 


0.22809 


0.24247 


0.25600 


0.18816 


c 




1.67676 


1.63743 


1.11252 


1.05000 


0.99691 


0.81641 


A = 


2 


3.72565 


4.06704 


4.61899 


5.43127 


6.30769 


7.22683 


m z 







1 


1 


2 


2 


3 


B 




0.49791 


0.29039 


0.30813 


0.24041 


0.25184 


0.20584 


c 




2.21655 


1.64297 


1.55874 


1.26233 


1.20167 


1.04900 


A = 


3 


4.32576 


4.60635 


5.23732 


6.01263 


6.89004 


7.81385 


m z 







1 


1 


2 


3 


4 


B 




0.54243 


0.34137 


0.36669 


0.28354 


0.23794 


0.21288 


c 




2.57492 


1.98559 


1.88097 


1.54922 


1.34215 


1.18987 


A = 


4 


4.85637 


5.11233 


5.73655 


6.53525 


7.41601 


8.33875 


m z 







1 


2 


3 


4 


5 


B 




0.56603 


0.37934 


0.30103 


0.25770 


0.22591 


0.20806 


c 




2.85052 


2.26315 


1.88523 


1.62705 


1.45132 


1.30691 


A = 


5 


5.34141 


5.59088 


6.21518 


7.01722 


7.90111 


8.82282 


m z 







1 


2 


3 


4 


6 


B 




0.57976 


0.40862 


0.33005 


0.28447 


0.25575 


0.20855 


c 




3.07820 


2.49976 


2.09758 


1.81530 


1.60855 


1.38742 


A = 


6 


5.79354 


6.04652 


6.68025 


7.46785 


8.34532 


9.27058 


m z 







1 


2 


4 


5 


6 


B 




0.58818 


0.43192 


0.35512 


0.26257 


0.24321 


0.22777 


c 




3.27473 


2.70721 


2.28554 


1.85449 


1.66757 


1.52332 



other resulting in a reduced mean square distance. However, in presence of Coulomb correlations there are values of 
A and 7 where electrons find energetically favorable to jump to outer orbits ( increasing the angular momentum and 
their mean square distance, as well) despite the effect of the magnetic field. A manifestation of this behavior comes 
in the form of jumps of a 2 (\pi — P2\ 2 ) for magnetic field 7=1 relative to the zero magnetic field case (7 = 0) as 
seen in Fig. [T] For instance for Coulomb correlation values A = and A = 1 electrons are squized closer to each other 
in presence of a magnetic field ( the values of a 2 (\p\ — P2\ 2 ) at 7 = 1 represented by empty circles are below the 
corresponding values at 7 = represented by filled circles), however for a larger correlation, such as A = 2 this is not 
the case any more. 

The increase in the mean square distance between the two electrons reflects the formation of a rotating clec- 
tron(Wigncr) molecule (REM or RWM). Quantum mechanically one can immediately see that (R) — 0, where R is 
the center of mass (CM) position and (pi) = —(p2)- This represents the quantum counterpart of the classical minu- 
mum energy configuration for two point charges in a harmonic trap that requires pi = —f>2- A plot of the electronic 
density will reveal that the electrons are mainly localized on a thin ring centered at the zero of the parabolic potential. 
The electronic density will exhibit the circular symmetry of the Hamiltonian and will only depend on distances (no 
angular dependence). The electronic density is insensitive to angular correlations which are very important, partic- 
ularly for large values of the angular momentum. A better probe of the angular characteristcs of the wave function 
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TABLE III: The optimal variational value of a 2 {\pi — fcf) for 2D quantum dot helium subject to a magnetic field as 



.2 II ~ |2\ 

a function of dimensionless Coulomb coupling parameter, A = e 2 a/(A-Keoe r hujo) = 0, 1, . . . 10 and values of magnetic field, 
7 = u c /ujo = 0,1,.. .5. The parameter a = WrmJo/Ti has the dimensionality of an inverse length. 





7 = 


7 = 1 


7 = 2 


7 = 3 


7 = 4 


7 = 5 


A = 


2.00000 


1.78885 


1.41421 


1.10940 


0.894427 


0.742781 


A = 1 


3.18193 


2.79350 


3.19963 


2.47708 


1.97594 


2.33444 


A = 2 


4.15428 


4.62431 


3.56731 


3.71459 


2.96353 


3.14825 


A = 3 


4.97151 


5.12920 


3.92383 


3.90825 


3.92870 


3.94601 


A = 4 


5.69236 


5.61781 


5.35004 


5.08275 


4.88226 


4.73529 


A = 5 


6.34774 


6.09048 


5.62186 


5.24316 


4.98332 


5.51752 


A = 6 


6.95629 


6.54753 


5.89000 


6.39368 


5.91907 


5.58307 


A = 7 


7.52912 


6.98994 


7.27192 


6.53384 


6.85216 


6.35681 


A = 8 


8.07359 


7.41855 


7.49712 


7.67346 


6.93629 


7.12900 


A = 9 


8.59467 


8.80371 


7.72156 


7.79918 


7.86399 


7.18439 


A = 10 


9.09607 


9.16053 


9.08988 


7.92518 


7.94189 


7.95454 
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FIG. 1: Variational mean square distance between two electrons, a 2 (\pi — p*2 1 2 ) for 2D quantum dot helium system in a 
perpendicular magnetic field as a function of dimensionless Coulomb coupling parameter, A = e 2 a/{^Tveoe r hu>o) for two values 
of magnetic field corresponding to 7 = u c /u)o = and 1. The line joining the data points at zero magnetic field serves as a 
guide to the eye. 



is the conditional probablity distribution (CPD) function [52| which for N electrons is quite generally defined as: 
P(p,po) — J2i=i $(pi ~ P) 3{Pj ~ where $ is the wave function under consideration. For the 

case of quantum dot helium (N = 2) this becomes: 



P(pJo) 



\*(p,Po)\ 
(*l*> 



(9) 



where the wave function is given from Eq. |(2} . When calculating the CPD function the position vector of one electron, 
pa is fixed, while p is moved so the resulting function of p measures the probability of finding one electron at p given 
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FIG. 2: Contour plots of the CPD function, P(p, po) corresponding to the displaced Gaussian ground state wave function for 
2D quantum dot helium at zero magnetic field (7 = 0) and A = 10. The black dot denotes the location of the fixed electron 
sittuated at position a po = (a 20 7^ 0,0). Distances are given in dimensionless units where a is the inverse oscillator length, 

yj niLUo/h. 



that there is one located at po. Obviously, for any choice of po 7^ the CPD function enables us to obtain the angular 
distribution of the second electron. Fig. [2] shows the CPD function for the ground state of 2D quantum dot helium at 
7 = and A = 10. The contour plots of the CPD function shown in the present work are calculated with p on the x 
axis with p corresponding to the distance at which the electronic density function (that is circularly ring shaped) has 
the maximum (crudely this can be thought as the radius of the thin ring). The black dot indicates pq. One can see 
clearly that the second electron is mainly localized in the opposite position of the fixed electron. This demonstrates 
that the quantum ground state has the same symmetry as the classical lowest energy configuration. 



DISCUSSION 



Beyond the ground-state energetics and ground-state angular momenta, we test the accuracy of the trial wave 
function in two well known limits: (i) infinite magnetic field limit (7 — > 00) where the ground state energy should 
approach the energy of a classical system of N = 2 point charges in a parabolic potential (adjusted by the quantum 
zero-point energy when Coulomb correlations are absent) and (ii) lowest Landau level limit where the ground state 
energy of the trial wave function must coincide with the lowest Landau level Laughlin wave function for N = 2 
electrons without the parabolic potential confinement. 

To study limit (i) we calculate the lowest classical energy, E c for N = 2 point charges in a harmonic potential which 
in dimensionless units is: 

^=7^t = |(2A) 2 / 3 . (10) 

We note that the classical energy is not afected by the presence or absence of a magnetic field. The classical ground 
state configuration for N = 2 electrons is one in which the respective positions of the particles are exactly opposite 
to each other at an optimal distance, p\ = —pi 7^ 0. 

Naturally, one cannot immediately compare the quantum variational energy, e to its classical counterpart, e c since 
without Coulomb correlations (A = 0) the lowest quantum energy, e is nonzero while the lowest classical energy, e c 
is zero. Such difference in energy between the two quantities represents the quantum zero-point energy (at A = 0) 
which in dimensionless units is: 
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FIG. 3: Variational ground state energy of 2D quantum dot helium in a perpendicular magnetic field, e — E/(huo) as a 
function of dimensionless Coulomb coupling parameter, A = e 2 a/ (4neoe r hLJo) for values of magnetic field corresponding to 
7 = ljc/uiq — 0,2,4, and 6. The solid lines represents the "adjusted" classical energy, e c + eo as a function of A calculated at 
the given 7 values. 



nujo V 4 

If we adjust the classical energy by eo the quantities to compare are e versus (e c + eo). Fig. [3] shows the variational 
ground state energy of 2D quantum dot helium in a perpendicular magnetic field, e and the adjusted classical energy, 
e c + £0 (solid lines) as a function of Coulomb correlation parameter, A for selected values of the magnetic field 
parameter, 7. Quite surprisingly there is very good agreement between the "quantum" ground state energy and the 
"adjusted classical" ground state energy at all magnetic fields including weak magnetic fields. As the magnetic field 
grows (increasing values of 7) the agreement only improves as can clearly be seen from the data. We also checked 
that in the infinite magnetic field limit (7 — > 00) both ^{variational) — eo and e(diag) — eo tend to e c (classical). This 
behavior was first noted by Yannouleas and Landman [22| in a study in which RWM wave functions were used to 
describe few electron quantum dots. They also recognized the importance of such finding in challenging the composite 
fermion picture of quantum dots which instead implies that e — eo— » as 7 - > 00. The same result is obtained here 
using not the RWM wave function, but the displaced Gaussian variational wave function. 

To check limiting behaviour (ii) we use Laughlin wave function to calculate the energies for the same A-s and 
7-s in which the displaced Gaussian trial wave function was used and then compare the results. The ground state 
energies obtained with the Laughlin wave function are displayed in Table [TV] The angular momentum, m z for which 
the lowest Laughlin energy is obtained is also specified. As expected, one notes that Laughlin wave function is not 
a good description of the system at weak (and zero) magnetic fields (for instance, there are several occasions in 
which the Laughlin ground state energy has the wrong angular momentum, such as the cases 7 = and A = 2, 3, . . . 
etc.). However, comparing the energies in Table ITT1 to Table HVl one sees that, in the limit of strong magnetic fields 
(increasing 7-s), the ground state energy of the displaced Gaussian variational wave function quickly approaches the 
Laughlin values. This is clearly seen in Fig. [4] where we plot the ground state energy corresponding to the displaced 
Gaussian variational wave function and the Laughlin wave function for A = 2 and 6 as a function of the magnetic 
field parameter, 7. The convergence of variational and Laughlin ground state energies in the limit of strong magnetic 
field is quite general and happens at any arbitrary value of the Coulomb correlation strength, A, though for clarity of 
the plot in Fig. [4] we display only the curves corresponding to A = 2 and 6. 
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TABLE IV: Ground state energies, e = E/(ftu>o) corresponding to the Laughlin wave function for 2D quantum dot helium 
subject to a perpendicular magnetic field for given values of dimensionless parameters A and 7 The angular momentum, m z of 
the ground state is also specified. 





7 = 


7 = 1 


7 = 2 


7 = 3 


7 = 4 


7 = 5 


A = 


2.00000 

m z = 


2.23607 

m z = 


2.82843 
m z = 


3.60555 

m z = 


4.47214 

m z = 


5.38516 

m z = 


A = 1 


3.25331 

m z = 


3.51671 

m z = 1 


3.98787 
m z = 1 


4.74972 
m z = 1 


5.64527 
ra z = 1 


6.54155 

m z = 2 


A = 2 


4.25331 
m z = 1 


4.17932 

m z = 1 


4.73309 

m z = 1 


5.47320 

m z = 2 


6.34988 
m z = 2 


7.24827 

m z = 3 


A = 3 


4.87997 
m z = 1 


4.84193 
m z = 1 


5.33361 
m z — 2 


6.09150 

m z = 3 


6.93735 

m z = 3 


7.84253 

m z = 4 


A = 4 


5.50663 
m z = 1 


5.45996 

m z = 2 


5.89253 

m z = 2 


6.61737 

m z = 3 


7.46625 

m z = 4 


8.37252 

m z = 5 


A = 5 


6.13329 

m z = 1 


5.95692 

m z = 2 


6.39990 

m z = 3 


7.11735 

m z = 4 


7.95855 

m z = 5 


8.86033 

m z = 6 


A = 6 


6.75994 
m z = 1 


6.45388 

m z = 2 


6.86566 

m z = 3 


7.57749 

m z = 4 


8.41976 

m z = 5 


9.31802 

m z =7 




2 1 , , , , , , , 

1 2 3 4 5 6 

y 

FIG. 4: The ground state energy curves, e = E/ (fouio) as a function of the magnetic field parameter, 7 for the case of displaced 
Gaussian variational wave function (Var) and Laughlin's wave function (Laughlin) for two values of the Coulomb correlation 
strength: A = 2 and 6. The solid lines joining the data points serve as a guide to the eye. 

Given the very good performance of the trial wave function on describing 2D quantum dot helium in an arbitrary 
perpendicular magnetic field, we plan to generalize the treatment to larger quantum dots. In this case the liquid or 
solid (crystalline) character of larger quantum dots in a magnetic field is rather non trivial and crucially depends 
on the values of 7, A as well as density of the system. To that effect we would describe any 2D quantum dot in a 
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perpendicular magnetic field with a generalized trial wave function of the form: 

JV 



*n = II V(P^) x (* - z j) n l D^D^) X (S) , (12) 



i<3 



where N is the number of electrons in the dot, x(S') = x( s ii s 2j • • • sjv) is the spin function for iVj spin-up and 
ATj_ = A" — A^f spin-down electrons and the space wave function has a Jastrow-Slater form. The determinants Dt (&) 



and Di($) are Slater determinants for spin-up and spin-down electrons built out of: (i) FD orbitals [53j, [54J or 
(ii) Gaussian localized orbitals [55j (for a crystalline state, only). The displaced Gaussian pair correlation factor, 
J(pij) as specified in Eq.Q gurantees the quality of the wave function at all magnetic fields ranging from weak 
(and zero) to strong and the integer quantum number, n = 0, 1, 2, . . . which takes even/odd values for respectively 
antisymmetric/symmetric spin functions determines the overall parity of the space wave function as required by the 
Pauli's principle. In such a case, a full VMC simulation would be the method of choice and the optimized trial wave 
function can be further used as a guiding function for DMC calculations. 



CONCLUSIONS 



To conclude, we have introduced a very accurate trial wave function for 2D quantum dot helium in an arbitrary 
perpendicular magnetic field. A key novel clement of this description is a Jastrow pair correlation factor that has 
a displaced-Gaussian form and contains two variational parameters to optimize. The variational energies are in 
excellent agreement with exact numerical diagonalization calculations at any value of the perpendicular magnetic 
field including weak (and zero) or strong fields. In agreement with the RWM formalism we find that, for given values 
of the Coulomb correlation strength, the quantum ground state energy of 2D quantum dot helium at any value of the 
magnetic field is close to the value of the "adjusted" classical energy. Quantum and classical energies converge in the 
limit of infinite magnetic field. For a given Coulomb correlation strength and in the limit of infinite magnetic field, 
the energies of the displaced Gaussian trial wave function agree very well with the energies obtained from Laughlin's 
wave function, though we note that Laughlin's wave function is a poor description of the system for weak (zero) and 
intermediate magnetic fields. For weak (zero) and intermediate magnetic fields a Jastrow pair correlation factor of the 
nature studied in this work should be included in the total wave function in conjuction with the Laughlin or RWM 
component, which are most effective in high magnetic fields. A Jastrow pair correlation factor, such as the displaced 
Gaussian factor introduced in this study, is essential to provide an accurate description of the system at any value 
of the magnetic field not limited to high magnetic fields only. Other Jastrow pair correlation factors such as those 
constructed through Pade approximations with several variational parameters [HI also provide quite an accurate 
description of the system and compare favorably with exact diagonalization results. However, the displaced Gaussian 
pair correlation factor is the most intuitive and the simplest physical choice that guarantees a consistent and excellent 
description of quantum dot helium system at all magnetic fields ranging from weak (zero) to infinity. A generalization 
of this trial wave function for A^-electron quantum dots in a perpendicular magnetic field is also discussed. 
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